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ABSTRACT 

The gap between ground and first excited state of tlie quantum-meclianical dou- 
ble well is calculated using the Renormalization Group equations to the second 
order in the derivative expansion, obtained within a class of proper time regula- 
tors. Agreement with the exact results is obtained both in the strong and weak 
coupling regime. 
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One of the most representative problems in quantum mechanics is the tun- 
nelhng of a particle through a potential barrier. It does not have a classical 
counterpart and cannot be handled in the usual perturbative approach. The dou- 
ble well potential provides a typical example of quantum tunnelling and, when 
employed in field theory, it also represents the most elementary toy model for 
the spontaneous symmetry breaking phenomenon and the related phase transi- 
tions and therefore it is extremely relevant in particle physics and in cosmology. 
The spectrum of the quantum mechanical double well can of course be evaluated 
numerically and it is known that the energy gap AE between the first excited 
state and the ground state, which appears as a consequence of the tunnelling be- 
tween the two classical vacua, has a singular behavior in the quartic coupling A: 
/\E oc exp(— 1/A) [|T|. The exponential behavior is well reproduced by the dilute 

istanton gas calculation, /S.E = 2\j2\/2/ (7rA)exp(— l/(3v^A)), which however is 
reliable only for very small values of A and becomes soon very different from the 
exact value of /S.E as A grows. 

Some time ago it has been suggested 0, Q that /S.E could be evaluated by 
making use of the Renormalization Group (RG). The original formulation of the 
non-perturbative RG obtained through a recursion procedure is due to Wilson 
Q and since then a differential formulation of the RG flow equations has been 
developed The RG equations involve the full action of the problem con- 



sidered and determine its flow as a function of a momentum scale starting from 
the bare action defined at a high momentum scale k = A down to the infrared 
region ~ 0. In the — > limit the running action contains the effects of all the 
fluctuations below the scale A and can be identified with the effective action. 

In 0, ^ the equation for the full action was approximated by an equation 
for the potential, which corresponds to the lowest order approximation in a sys- 
tematic derivative expansion of the action. In a further approximation was 
considered by expanding the potential in a polynomial series and then truncating 
this series. The flow equation for the potential was then reduced to a set of cou- 
pled ordinary differential equations for the various coefficients of the series. In 
the full partial differential equation for the potential was instead studied with no 
further approximation. In both cases the results are practically the same: when k 
is lowered the potential evolves from the bare double well shape, reducing the bar- 
rier until it totally disappears. When k is close to zero the potential has become 
convex with the central region between the original minima almost fiat. This is 
already an interesting feature because the convexity of the effective potential is 



an exact property |TT| which can be observed by making use of non-perturbative 
approaches (see e.g. [|T^) and cannot be recovered in perturbation theory. On a 
more quantitative level, however the flow equation for the potential is not fully 
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satisfactory. In fact when it is used to compute the quantity AE, which as ex- 
plained in is related to the second derivative of the effective potential at the 
origin, it provides accurate results for large values of the coupling A, but when 
approaching A = the correct exponential behavior is not recovered. Therefore 
the regime in which the quantum tunnelling effects are relevant (A ~ 0) is not 
well described neither by the approximation considered in nor by the one in 

i- 

In this Letter we examine the natural improvement on the lowest order ap- 
proximation of the flow equations in the derivative expansion, and specifically we 
include the next term in this expansion, i.e. the coefficient of the kinetic term 
of the action, Z{k,x), which depends both on the running scale k and on the 
spatial coordinate x. In field theory this is nothing else than a wave function 
renormalization and in quantum mechanics it has the role of a position depen- 
dent effective mass. We expect that the inclusion of this term, which represents 
the first correction to the fully local information carried by the potential could 
be helpful in describing the complex dynamics of the tunnelling. 

Our aim is not only the analysis of a possible improvement in determining AE, 
due to the inclusion of Z{k,x). In fact we are particularly interested in checking 
the reliability of a particular version of the RG flow obtained by means of a 
regulator introduced in the Schwinger Proper Time (PT) formalism [l^, |15 



W\. As discussed in []TB[ and more in detail in [|T^] this particular flow, unlike the 
so called Exact Renormalization Group (ERG) flow, does not have a first principle 
derivation and there is no proof of convergence to the effective action in the 



A; — »• limit. However, as shown in [jl5|, |T6[, the PT RG equations are particularly 
interesting because they provide excellent determinations of the critical exponents 
at the non-gaussian fixed point in three dimensions, certainly comparable to the 
ones obtained with the ERG equations. The problem of determining AE in the 
quantum mechanical double well is therefore another good check for the PTRG 
because at least in this case the exact results are available. 



In [^, |T6| the flow equations of the PTRG are derived in the required approx- 
imation of the derivative expansion of the action, i.e. the coupled partial differen- 
tial flow equations for the potential and for the wave function renormalization are 
displayed for the particular PT regulator fk^sZk"^) = exp{—sZk^) J2iLo {sZk^y/i\ 
where s is the proper time and m is an integer index. Actually the best values for 



the critical exponents in |15, 16 1 are obtained in the limit m — > oo and therefore 



here we are particularly interested in this limit. In |]T6[ the limit m ^ cxd is for- 
mally taken in the fixed point equations by redefining the field variables. It is also 
noticed that the second derivative of the potential determined by the m-dependent 
equations behaves like 1/m for large m. This behavior would imply in the present 
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case that AE = in the large m hmit. However this problem can be avoided if 
one, instead of considering the field redefinition as in [jl6[, would simply start with 
a slightly different cut-off function fklsZk"^) = e:!q){—sZ'mk'^)Y^'^Q{sZmk'^Y/i\ 
where fc^ has been replaced by mk"^ . This is practically equivalent to what has 
been done in Jl^ where the same rescaling of k has been performed directly in 
the flow equations. However the latter approach mixes the two limits m ^ oo 
and — and this could generate some confusion. 

If the new cut-off function fk is used instead of /fc, a new set of m-dependent 
flow equations are obtained whose large m limit (this time with no redefinition of 
the field or of the scale k) yields again the equations which have been determined 



m 



IHl and 



The m-dependent flow equations obtained with the cut-off fk 
as we shall see below provide for the double well problem, in the limit m — oo 
a non- vanishing second derivative of the potential and therefore a finite value of 
AE. The procedure outlined in |[T^ to derive the flow equations, with the cut-off 
fk replaced by /fc, yields the two coupled partial differential equations for the 
running potential V{k, x) and the wave function renormalization x) 



dV 



Zk^ + V"/m 



dk ^ ' \Zk^ + V"/m 



m+l-D/2 



(m + 1 - D/2) 
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where each prime indicates a derivative w.r.t. the spatial coordinate x and the 
constant a is expressed in terms of gamma functions 

T{m + 1 - D/2) 



Eqs. 
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(47r)^/2r(m 
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Note that all equations are expressed in terms of dimensionful quantities (with 
^ = 1). The determination of AE requires the numerical integration of the 
coupled equations for V and Z with D = 1. We have to fix the initial conditions 
of the differential equations which we take at a value k = A, with A much larger 
than any other scale in the problem. The double well potential 

= ^M^x^ + Xx' (6) 

corresponds to the initial condition of the flow of the potential V{k, x) whereas, 
according to the normalization of the kinetic term in the bare action, we take 
Z = 1 as initial condition for Z{k, x). The scale of all the dimensionful quantities 
is fixed by the choice = — 1. We shall also consider, for comparison, the convex 
anharmonic oscillator, and in this case we choose = 1. Correspondingly the 
scale at which the initial conditions of the parameters are fixed is chosen A = 
1500. The numerical resolution of the two coupled partial differential equations 
is performed with the help of the NAG routines. 

In Fig. 1 V"{0, 0) and Z{0, 0) obtained from Eqs. (^0) are reported versus 
1/m. Namely black circles and diamonds correspond respectively to Z{0,0) and 
V^(0, 0) for = 1 and A = 0.4. White circles and diamonds correspond to 
Z{0,0) and V{0,0) for = —1 and A = 0.05. In all cases we have inserted at 
1/m = the corresponding values obtained from Eqs. (^,|^). The convergence to 
these latter values is clear. Therefore, as we had anticipated, in the limit m — cx), 
Eqs.(|I|J^) provide non- vanishing values of the second derivative of the potential 
which converge to those obtained from Eqs.(^J^). 

Since we are particularly interested in the 1/m = case, from now on we 
consider Eqs.(0J^). The qualitative behavior of the running potential along the 
flow is explained in detail in 0] and we observe a similar trend. In Fig. 2 
the second derivative of the running potential (continuous lines) and the wave 
function renormalization (dashed lines) obtained at various values of k (namely 
k = 1500,10,1,0.5,0.2,0.1,0) for = -1 and A = 0.06 are plotted. Below 
k = 0.2 the second derivative of the potential at the origin becomes positive and 
the potential convex. At the same time Z, which for most of the running stays 
close to one, rapidly increases in the region between the classical minima with a 
sharp fall to one outside that region. The highest peak corresponds to the final 
output at k = 0. 

Finally we turn to the evaluation of the energy gap AE. As explained in 
as long as the potential is analysed keeping Z = 1 fixed during the flow, AE is 
to be identified with ^V"{0, 0), but in the approximation here considered with a 
running wave function renormalization, the energy gap AE corresponds, in field 
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theory language, to the renormahzed mass, i.e. 



AE = ^V"iO,0)/Z{0,0) (7) 

Therefore we determined the values of AE from Eqs.(^,^) according to Eq. (|^) 
(these results are indicated as AEpTnio) and, for comparison, we also evaluated 
AE to the lowest order in the derivative expansion (indicated as AEpTio), i-e. 
solving the flow in Eqs.(^,^ keeping the wave function renormalization constant 
Z = 1. The energy gap to the lowest order in the derivative expansion with 
the flow equation used in 0, which is identical to the Wegner- Houghton local 
potential equation in D = 1 (see ^) 

has also been calculated (AEy^h) and all these results are eventually compared 
with the exact estimate AEexact, derived by numerically solving the quantum 
mechanical eigenvalue problem. In Table 1 these estimates of AE are collected 
for various values of the coupling A both for = 1 and = — 1. The values 
found for Z{0, 0) are also included. From the accuracy of the numerical code the 
estimated errors in Table 1 are at most two or three units on the last digit. 

In the = 1 expected from 0, |^, everything is rather smooth. In 

fact we observe that the differences among the various determinations of AE are 
extremely small and Z{0,0) is practically equal to its bare value Z{A,x) = 1. 
However little discrepancies appear at large values of A (strong coupling regime) 
and AEyjh turns out to be a slightly better approximation than AEpnpio and to 
be comparable to AEpTnio- 

The double well problem with = — 1 is rather different. Here, from p|, Q we 
know that the RG equations work better in the strong coupling regime. Actually 
this feature is true for AE^h which, in the region where the tunnelling plays 
an important role (A < 0.1), does not reproduce the exponential behavior of 
AEexact- On the other side, AEppio which is less accurate at A = 0.4, turns out 
to be more reliable at small values of A. But the remarkable result is that the 
higher order approximation is very accurate. In fact the correction due to the 
wave function renormalization increases when A is decreased and the agreement 
between AEpTnio and AE^xact is always excellent. Unfortunately below A = 0.05 
the peak in Z(0, x) shown in Fig. 2 becomes very narrow, generating problems in 
the numerical resolution of the flow equations and we were not able to determine 
the corresponding values of Z(0, 0) and of the energy gap. In any case for such 
small values of A the numerical errors are no longer negligible and a more refined 
analysis should be performed. Even an approximated analytical investigation of 
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this regime, similar to the one discussed in 0, would be useful, although the 
problem of two coupled equations is much more difficult to handle. 

In conclusion we have determined the energy gap between the ground state 
and the first excited state of the quantum mechanical potential Vdw in Eq- (|§) 
both for positive (quartic anharmonic potential) and for negative (double well 
potential) curvature M^, by making use of the RG techniques. In particular we 
used a Proper Time regulated version of the RG flow equations in the derivative 
expansion approximation, truncated to the second order. The results show the 
reliability of the PT RG flow which, even to the lowest order provides good 
estimates of the energy gap and, above all, they show the sensible improvement 
coming from the inclusion of the wave function renormalization in the quantitative 
analysis of a tunnelling process. 
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TABLE 1 



A 




^EpTlo 


^Egxact 


^EpTnlo 


Z(0,0) 


M2 = 1 


1.0 


1.9291 


1.9464 


1.9341 


1.9380 


1.0052 


0.4 


1.5450 


1.5556 


1.5482 


1.5498 


1.0037 


0.1 


1.2091 


1.2127 


1.2104 


1.2109 


1.0013 


0.05 


1.1201 


1.1218 


1.1208 


1.1210 


1.0006 


0.03 


1.0774 


1.0784 


1.0779 


1.0780 


1.0003 


0.02 


1.0538 


1.0544 


1.0540 


1.0542 


1.0002 


M2 = -1 


0.4 


0.9654 


0.9897 


0.9667 


0.9730 


1.0217 


0.3 


0.8173 


0.8404 


0.8166 


0.8233 


1.0273 


0.2 


0.6212 


0.6416 


0.6159 


0.6227 


1.0416 


0.1 


0.3297 


0.3280 


0.2969 


0.3027 


1.1321 


0.07 


0.2238 


0.1848 


0.1539 


0.1562 


1.3343 


0.06 


0.1902 


0.1311 


0.1031 


0.1028 


1.5548 


0.05 


0.1576 


0.0806 


0.0562 


0.0532 


2.1270 


0.04 


0.1259 


0.0496 


0.0210 






0.03 


0.0947 


0.0329 


0.0036 






0.02 


0.0637 


0.0201 


0.0003 







Tabic 1: Various determinations of AE, including its exact estimate, and of 
Z(0, 0) for some values of the coupling A, both in the = 1 and in the = —1 
case. 
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Figure 1: V"{0,0) ( black and white diamonds) and Z{0,0) (black and white 
circles) obtained from Eqs.(|l],^ respectively for = 1, A = 0.4 and for = 
— 1, A = 0.05, vs 1/m. The limiting values at 1/m = are obtained from Eqs. 



9 







1 



2 



X 

Figure 2: V"{k,x) (continuous lines) and Z{k,x) (dashed lines) plotted vs x 
at k = 1500,10,1,0.5,0.2,0.1,0, as obtained from Eqs. (|,|) for = -1 and 
A = 0.06. 
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